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In the framework of perturbative QCD, double inclusive cross sections for high pt 

■ parton production in nucleus-nucleus collisions are computed with multiple rescat- 

■ tering taken into account. The induced long-range correlations between numbers of 
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jets at forward and backward rapidities are found to reach 10 20% for light nuclei 
at ^/s = 200 GeV/c and to be suppressed for heavy nuclei and LHC energies. 

I. INTRODUCTION 

Particle production in high-energy heavy ion collisions is now at the center of exper- 
imental efforts to discover the quark-gluon plasma The observed particle spectra 



^ I are the result of different mechanisms which are responsible for the creation of initial 

\^ • high-p( partons, their propagation and subsequent hadronization. (Experimentally the 

O . initially produced partons can be traced as jets of hadrons, so that in the following 

^ ! we shall often use the term 'jet' for the produced parton in the pure theoretical sense, 

neglecting all the subtleties related to its actual determination in the experiment.) To 
p • be able to see the formation of the quark-gluon plasma against the background of more 

conventional effects, such as as gluon emission due to bremsstrahlung and multiple hard 
Qj I collisions in the nuclear medium, one has to fully study the consequences of the latter. 

'. Jet_quenching due to gluon emission has been studied in considerable detail (see e.g. 

, — j)- Also the effects of rescattering in single particle inclusive spectra initiated by 
^ . the observation of [s*] has been considered in a series of articles for multiple soft 0] 

and multiple hard pQCD scatterings U, 13, 13 ■ This paper generalizes the 



study of multiple hard collisions to double inclusive cross-sections and the following 
long-range correlations between the secondaries. Observation of such correlations has 
always served a very precise instrument for analyzing the dynamics of the interaction. 

Following the framework introduced in we consider only relatively hard colli- 
sions which allow for the perturbative QCD approach. The colliding nuclei are assumed 
to contain a variable number of partons with initially small transverse momenta, which 
become large as a result of hard collisions between partons belonging to the projec- 
tile and target. The latter are assumed to move fast along the collision axis in the 
opposite directions. In the present work we restrict ourselves to the study of high 
Pt parton production at rapidities well separated from central rapidity region. Here 
longitudinal momenta of partons in the final state are much larger than the transverse 
ones. So we neglect attenuation of the former during the collisions and assume them to 
be conserved throughout the nucleus-nucleus interaction in agreement with the stan- 
dard Glauber treatment. We neglect intrinsic correlations between partons inside the 
colliding nuclei which correspond to nuclear shadowing and take the partonic distri- 
butions just as a product of such distributions for the nucleons smeared out with the 
standard nuclear profile functions. This implies that we correctly describe quantum 
evolution of these distributions except for effects coming from interaction of partons 
between different nucleons in the nucleus, which should be studied from the DGLAP 
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evolution equation for the nucleus as a whole or, at very small x, in the framework 
of the BFKL evolution (e.g. in the Colour Glass Condensate approach, see and 
references therein). Jet quenching and hard parton contents of the participants are 
some of other important phenomena to be taken into account in the full treatment of 
jet production. Thus the physical phenomenon we study is restricted to multiple hard 
elastic scattering between partons. Our study is then to be considered as a baseline 
calculation to be followed by inclusion of all the above-mentioned effects 



Our formalism closely follows to which paper we refer the reader for details. 
As in our basic approximations in the study of multiple partonic collisions are 

1) purely elastic partonic collisions with conservation of longitudinal momenta and 

2) factorization of the S'-matrix into the product of elementary partonic S- matrices 
(the Glauber approximation). We shall also use the notations of To make our 
presentation more self-contained we reproduce some basic notations below. 

The initial states of the colliding nuclei A and B are represented as a superposition 
of states with different number of partons, having specific values of scaling variables 
and transverse positions with respect to the target nucleus center. We combine them 
into a single argument z = {v, b} for nucleus A and u = {w, c} for nucleus B. The 
(variable) number of partons will be denoted correspondingly by n and /. In these 
notations the initial state of the nucleus A is 

oo „ 

\A) = ^ / dTA{n)'^A,n{^i^- ■ ■ ^^n)\n,Zi), (1) 



n=l 



where dTA{n) = Y[ d^Zi stands for the phase space volume of the configuration. Wave 

1=1 

functions A,n of the n-parton configuration, symmetric in their arguments, are nor- 
malized to fulfill {A\A) = 1: 



y"rfr(n)|^A,„(^i,...,^„)l% J]W^n = l. (2) 



We will assume a Poissonian distribution for the number of partons: Wn = e~^"'*/n! 
and, as we neglect intrinsic partonic correlations, a factorization property for the wave 
function 

n 

A,n{Zi)\^ = CnWpA{Zi). (3) 

i=l 

The Poissonian distribution immediately implies 

Cn = ^e-<">, (n) = j d'zpAiz), (4) 

where (n) corresponds to the average number of partons in nucleus A. The same nota- 
tions are used also for nucleus B. 

The paper is organized as follows. In sections |n] and IIIII we derive expressions for 
2-parton inclusive cross section for two different cases. The case when both partons 
belong to the same forward rapidity region is examined in the section m derivation 
for the forward-backward case is presented in section IIIII The form of the two-parton 
cross sections implies that non-trivial forward-backward correlations emerge in high-p^ 
parton production. In section IIVI a quantity measuring the strength of correlations. 
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correlation coefficient, is proposed which appears to be expressed in a simple way in 
terms of the derived cross sections. Numerical values for the correlation coefficient 
followed by discussion are presented in section |3for light and heavy nuclei interactions 
at RHIC and LHC energies. The interpretation of our results and conclusions follow 
in section IVII 



II. TWO JETS IN THE SAME DIRECTION 



The double inclusive cross-section to produce two jets in the same (forward) direction 
is obtained as a direct generalization of the single inclusive cross-section in |11|]. As 
indicated above, here we limit our studies to the production of partons at rapidities well 
separated from zero and neglect changes in the parton longitudinal momenta during the 
interaction. This implies that both observed partons originate from the same nucleus 
(projectile A). 

We have to fix the final state A' of the nucleus A to describe the two produced 
partons. Since the total wave function is symmetric in all the partons we choose 
th em to be the ffist and the second one and account for other possibilities by factor 
y/n{n - 1): 



-^A',n{^-i-^^2,Z3,...Zn) = n{n - 1) ^^^^ (zi ) ^^^^ (^2) ^A',n-2 (^3, • • -^^n) 



(5) 



where a' ,n-2{zz-, ■■■Zn) is the symmetrized wave function of the unobserved n — 2 
partons and tp^iizi), i = 1,2 are the wave functions of the observed partons in the final 
state, ai combine scaling variables and transverse momenta pi and p2 of the latter. 

The probability to observe the two partons is given by the modulus squared of the 
S'-matrix element {AB\S\A'B') summed over all possible final states of the unobserved 
partons. In our approximation (purely elastic collisions) the number of partons is not 
changed by the interaction. So the S'-matrix is diagonal in the basis {n, zf, I, Ui}: 

{n',z[,l',u'AS\n, Zi;l,Uj) 



Snn'Su'll5^'\z, - z'^l[5^'\ui - u'^SUzl, • • ■ , z^\u,, . . . , Ui) . (6) 



i=l 



Following ref. we take the square modulus of (6) making use of the specific form 
for the final-state wave function of the projectile nucleus (0). The double differential 
cross-section at a fixed overall impact parameter jS which follows reads 



da, 



a\a2 



d^(3 



oo „ 

^^ra(ra-l) / dzidz2dz[dz2iJaAz[)iJa2{z[)'ilj*^^{zi)'ipl^{z2) 



n=2 l>0 



dTA{n - 2)rfrij(/)^^„(zi,4,2;3, ...Zn)'^ A,n{zi, Z2, 2:3, ...2:„)|^b,/(mi, • • • ,^i«)P 



Snl{Zl-, 2^2, Z'i...Zn\Ui, ...Ui) 



S*nl{z[, Z'^, Z:i...Zn\Ui, ...Ui) 



(7) 



The product of square brackets in ((Tj) gives four terms. The double inclusive cross- 
section for transverse momenta pi,P2 3> 1/Ra,b corresponds only to the term which 
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originates from the product of the S'-matrices. Indeed in the Glauber approximation 
we assume that the S'-matrix is a product of S'-matrices for pair parton colhsions. Then 
terms hnear in S or S* diagrammatically correspond to cutting the forward scattering 
amphtude either to the extreme left or to the extreme right of all partonic interactions. 
In both cases the momenta of intermediate partons in the cut coincide with their values 
in the initial colliding nuclei, so that their transverse momenta are small, of the order 
of typical nuclear scale 1/Ra,b (see for more detail). 

We present each of the S'-matrices as a product of the elementary partonic ones, 

n I 

Srd{Zl, . . . , Zn\Ui, ,Ui) = YIYI 

Here Sij = 1 + ia{zi,Uj) and a{zi,Uj) are the parton-parton scattering matrix and 
amplitude respectively. Since only elastic parton scattering are considered, unitarity 
of the partonic s-matrices means SijS*j = 1. So we obtain that 

^nli^l^ ^25 Zs...Zn\Ui, ...Ui)Snl{Zi, Z2, Z3...Zn\Ui, ...Ui) 



I 

= Y[Y[[^ + H<,u,)]*[l + ta{zi,uj)]. (9) 

i=l,2j=l 

and does not depend on z^, ...Zn- This allows to integrate over these variables and sum 
over n to produce the density matrix of the nucleus A for a pair of partons: 



'^n{n - 1) / dTA{n - 2)^!\^^{z[, z'^, z^, ...2;„)^A,n(2i, 2:2, ^3, ■■■Zn) = Pa(^i, ^2^1, 4) 

n 

(10) 

Factor n(n — 1) again accounts the possibility that any two of the partons entering the 
symmetric wave function can be chosen. Thus we get 

—^^ = ^ j rf2;irfz2rf2;id4^/'„,(zi)^/'„2(zi)^;^(zi)^;2(^2)PA(2;i,2;2ki,4) 



dTB{l)\^BAuj)?\ n \{[^ + ici{zluj)]*[l + ia{zi,Uj)]-l\. 
Making use of the factorization of the wave function Q we have 



(11) 



2 1 - ' 

Pa{zi,Z2\z[,Z2) = Pa{zi\z[)pa{z2\z2) and \^ sA'^j)? = jf ^^^WPB{uj\uj). (12) 

i=i 

Assuming factorization of partonic transverse and longitudinal degrees of freedom, for 
equal scaling variables Vi = v[ = have 

PA{vM%)=PA{v)pA{h%), (13) 

where Pa{v) is the mean parton number distribution and pA{hi\h[) is the transverse 
part of single parton density matrix. For equal arguments it goes into the standard 
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profile function of the projectile nucleus PA{bi\bi) —TA{bi — P) (recall that the origin in 
the transverse plane is in the center of the target nucleus B). Similarly for the target 
nucleus pB{uj\uj) = PB(wj)TB{cj). 

As a result of this factorization the inclusive cross section transforms into 



lAA{p,yiPi,y2,P2) 



dyi dy2 cPpi d?P2 d'^P 



yi ,2/2 >o 



/ n (^d%,d%[e^^^('^-''^^PAiv,)pAm) 



1=1,2 

I 



I \{d^c^dw^PB{w^)TB{cA n \{[l + ia{z[,u^)]*[l + ia{z,,u^)]-l\ (14) 

where the rapidities yi and y2 correspond to the scaling variables Vi and V2. The 
subscript AA for the inclusive cross-section indicates that both jets are produced from 
the nucleus A, although the collision is between A and B. 

The non-trivial integral over Uj — {wj, Cj} factorizes to give the Z-th power of 

J{vi,V2,bi,b[,b2,b2) ^ / d^cdwTB{c)PB{w) JJ^ [l+ia(i;i, 6--c)]*[l+ia(i;i, fej-c)] 

i=l,2 

(15) 

so that after summation over I we get 

/AA(/3,yiPi,y2,P2) = j n (rf'&irf'&:e^^'(''-''^)p^(6,|69^(^e-^(-'-'''^'''-''^'^^ 

(16) 

We are left with a problem of calculation J, which is now more complicated than for 
the single inclusive cross-section. Prom the 16 terms in the product in (7) we separate 
a part 



)-{i) _ 1 



i+E(i 

1=1,2 ^ 



ia{yi, w; b[ — c)]* -|- ia{vi, w; bi — c) + [ia{vi, w; b[ — c)]*ia{vi, w; bi — c) 



Physically it corresponds to the case when either their is no partonic interaction at 
all or only one parton from nucleus A interacts with a given parton from nucleus B. 
Integration over c and w of this part repeats that for the single inclusive case and gives 



where 



Fb{x, b) = I dwPBiw) J d'p^^^^e^^' (18) 

Putting (9) into (8) we find that the corresponding part of the double inclusive cross- 
section factorizes into a product of two single inclusive ones: 

lAAi(^^yi^Pl^y2,P2) = lA{yi,Pl)lA{y2,P2) (19) 
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where lA{y,p) is the inclusive cross-section to produce a parton from nucleus A in its 
collision with nucleus B at impact parameter (3 (the latter dependence implicit): 

I^^y^p) = p^fff = P^(v) [ d%d'rTAib~(3)e'^^( e^BW(F,(.,r)-F,(.,o))_g~T,WF,(.,o)^ 
dyd^pd^P J \ J 

(20) 

So the part corresponds to independent production of the two partons as expected. 

The rest 9 terms correspond to the case when both the observed partons interact 
with the same parton from nucleus B. Evidently in this case the two observed partons 
have to be located close to each other, since the partonic interactions are short-ranged. 
As a result we shall find under a single integral over the impact parameter b a square 
of the profile function T'^{b). Thus this part of the inclusive cross-section will have 
the order ~ \j A^l'^ as compared to the independent production part. Having in mind 
calculation of correlations, we can neglect this part in the first approximation, as we 
shall see in the following. 



III. TWO JETS IN OPPOSITE DIRECTIONS 

In this case the wave functions of the observed partons V'ai,a2 from nucleus A (moving 
along the z-axis) and nucleus B (moving in opposite direction) will be denoted by 
'^^{z) and i^qiu) respectively, where p and q combine their scaling variables [v and w 
respectively) and transverse momenta. The n-particle wave function for the final state 
of the nucleus A takes the form: 

and similarly for the nucleus B. 

Again, writing down the relevant matrix element for the transition amplitude and 
again following ^J], we obtain the double inclusive cross section at a fixed overall 
impact parameter /? 

lAB{P,yi,p,y2,q) = ^^^'^^'^ / dT'^{n')dT'B{l')dTA{n)dTB{l) 

nln'l'>l A'B' 

^Xn-l(^-2---^n)^*B,l-li'>^2---Ul)'^AA^l---Zn)^BA'^l---Ul) 

[Snli4---ZnWvU'i) - l][Snl{Zi...Zn\Ui...Ui) - 1] (21) 

where the set Zi = {vi, bi} of longitudinal momenta and transverse positions corre- 
sponds to the nucleus A, whereas the set of Uj = {wj, Cj} to nucleus B. The subscript 
AB indicates that the two jets are now produced from different nuclei. 

Summing over all unobserved states of final nuclei A' and B', we can perform inte- 
grations over Z2...z'^ and u^.-u'.^ to obtain 

lAB{P,yi,P,y2, q) = / dzidz[dTA{n — l)duidu[drB{l — 1) 

nl>l 



^A,n(^l'^2, ...Zn)'^*B^iiu[,U2...Ui)[S*{z[,Z2...Zn,u[,U2, ...Ui) - l]^pp{z[)ipq{u[) 
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^p*{zi)tp*{ui)[S{zi...Zn,Ui...Ui) - 1]"^ aA^I-'-Zu)"^ B,l{^l--Ul) (22) 

The product of the amphtude and its conjugate gives three terms 

[S* -1][S-1] = [S*S - 1] - [S* -1]-[S-1] (23) 

For the same reason as for the case of two jets in the same direction discussed in the 
previous section, the second and third terms in ()23j) give the partonic spectrum only 
at very small values of transverse momenta and can be neglected. To be more explicit, 
consider terms with no more than one partonic collision with the projectile and/or 
target. Inserting transverse parts of the wave functions of observed particles in the 
form ip ~ e^^^ we arrive at the integral (for the term with S) 

J d%i£cd^CidVc[d\e'P^'+''^^'-'P^'^-'<X{b, h, c, ci) 

*^(6i..6..)vI/b(ci..c..)v1/;:i(6;..6..)v1/;^(c;..c..) (24) 

where 

X = ia{bi — c)ia{ci — b) , or ia{bi—c), or ia{b — ci) (25) 

The scale at which the wave functions of nuclei change is much greater than 1/p or 
1/g so all of these terms after integrations in 6'i,c'^ give contributions proportional to 
the product of delta-functions S{p)6{q). If we are interested in production of particles 
with large transverse momenta all these contributions can be neglected, so that in 
only the first term contributes. 

We present the total S'-matrix elements as 

/ n n,l 

Snl{zi..Zn\Ui..Ui) = JJ Sij JJ Sn JJ 

j=l 1=2 k,m=2 

Due to unitarity of the partonic s-matrix, SijS*j = 1, the product S*S takes the form 
S^i{z[, Z2..Zn\u[,U2..ui)Sni{zi, Z2..Zn\ui, U2..ui) = (1 + ia{z[,u\)y{l + ia{zi,ui)) 



I n I n 

JJ(1 + ia{z[,Uj)y JJ(1 + ia{zi, u\)y JJ(1 + ia{zi,Uj)) JJ(1 + ia{z,, m)) (26) 

j=2 i=2 j=2 1=2 



To move further, in contrast to the derivation in we are forced from the start 
to make the assumption of factorization of the nuclear wave functions. 

^^(zj, Z2...Zn)'^A{zi, Z2...Zn) = p a{z[\Zi) Pa{z2\z2) ■■ PA{Zn\Zn) (27) 

and similarly for nucleus B, where the density matrices pA{zi\z'^ have the factorization 
property (fT!?jl and normalization (jH). 
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Assuming for simplicity a central collision (in other case the argument in all T^s 
must be shifted by the impact parameter (3) we get 

lAB{f3,yi,p,y2,q) = ^^^ dzidz[duidu[tpp{z[)ijq{u[)ij*{zi)tlj*{ui)pA{zi\z[)pB{ui\u[) 



nl>l 



[1 + ia{z[, u[)]*[l + ia{zi,ui 



)]( n / dujPBiuj\uj)[l + ia{z[,Uj)]*[l + ia{zi,Uj)] 



JJ / dziPA{z.i\zi)[l + ia{zi,u[)]*[l +ia{zi,ui)] \ - / JJpbJJ 

1=2-^ ' 



Already at this stage it is convenient to separate from the total inclusive cross- 
section its part which does not contain interactions between the observed partons, that 
is the term which comes with unity from the product of the first two square brackets 
in (28). Comparing with |0| we find it to be: 

lAB{P^yi^P^y2,(l) = \ ^ I)] / dZidz[^p{z[)^*{zi)pA{Zi\z[) 

^ n>l ^ '' •' 



JJ^ / dujPB{uj\uj)[l + ia{z[,Uj)]*\l + ia{zi,Uj)]- / J^rf^jpsl 

3=2'' 3=2 

I 2^ ^l_iy_ J duidu['4j*{ui)'iljg{u[)pBiui\u[) 
JJ / dziPA{zi\zi)[l + ia{zi,u[)]*[l + ia{zi,Ui)]~ / Y[dZipA{zi\zi 

i=2 1=2 

= lA{yup)lB{y2.q) (29) 

where Ia{b) is the inclusive cross-section to produce a single jet from nucleus A{B) in 
AB collisions at impact parameter (3 (see (20)). So this term factorizes into a product of 
two independent single inclusive cross-sections to produce each of the observed partons. 
It is important that this result is exact and not based on the smallness of the parton 
interaction range on the nuclear scale. In particular, its validity is not spoiled by 
corrections of the order 1/ A^^^. 

Now we turn to the rest terms in (28). In them the integrals in round brackets 
can be done exactly as in the standard Glauber derivation, taking into account that 
the space range of the partonic interaction is much smaller than the nuclear scale (see 
jll|). After that summations in n and / can be easily performed and we get 

-^AB(/^5 2/i,P,l/2,g) = j d%d^b[d^cid^c[ 

e^^(''^-''^)+^''(^^-^'^W(t;, &i|&'i)pb(^, Clio 
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where 

Eb{v, h, b[,) = Tsiih + b[)/2)iFBiv, h - b[) - Fb{v, 0)) (31) 

Ea^ is defined by a similar formula for nucleus A and yi, y2 ,Pi, P2 are the rapidities 
and transverse momenta of the observed particles(partons), and v,w are longitudinal 
momentum fractions corresponding to yi and y2 respectively. Functions Fb{v, b) is the 
Fourier transform of the transverse momentum distributions I{v,w,p) in the elastic 
scattering of two partons with scaling variable v and w, averaged over the longitudinal 
momenta of nucleus B partons (see jlll|): 

FB{v,b) = I dwP{w) I ^^I{p,v,w)e'^' (32) 

It is convenient to introduce new variables for integration. Define 

n = b,- b\- r2 = c^- c\- b={bi + b\)/2- c={ci + c;)/2 
In these variables the cross section is rewritten as 

/S(/5,2/i,P,l/2,g) = j d\-,£r2dVce'P'^-^'''^'pA{v,b,\b[)pBiw,ci\c[) 



I (^[l+ia{b-c- ""-^-^^ V, w)Y[l+ia{b-c+'-^^, v, w)] - 1^ e^^("''''^^)+^^("''^''^^) -1 1 

(33) 

where in the new variables 

EBiv,b,n) =TB{b)iFBiv,n) - FBiv,0)) (34) 

and similarly for E^. 

Assuming that E{..b,c...) and density matrices change significantly only when b and 
c suffer macroscopic shifts of about nucleus radius, we can perform integrations in c. 
Take the terms in (20) containing single amplitudes ia{b — c ± (ri — r2)/2). On the 
nuclear scale they can be effectively substituted as 

ia{b - c ± (ri - r2)/2) ^ m{0)S^{b - c ± (n - r2)/2) ~ m{0)5^{b - c) 

where d is the amplitude in the transverse momentum space and we have used that 
(ri — r2)/2 is small on the nuclear scale. In the term with the product 

a{b-c+ in - r2)/2)a*{b - c - (n - r2)/2) 

we pass from the integration variable c to r = 6 — c + (ri — r2)/2 in which this product 
takes the form 

a{r)a*{r — ri + 
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Obviously r is also small on the nuclear scale, so that c ~ 6 and we may take all the 
c-dependent function out of the integral over r at c = 6. The integral over r takes the 
form 

/ d ra{r)a {r — ri + r2)— / d pe '^^^ ^' — — 

J J d^P 

In this way we finally get for the second part of the inclusive cross-section 



{(3,y,,p,y2,q)^PA(v)PB(w) j d%TA{h)TB{h) j d\^d\2e'^ 
(^jdne^^^^^-^^^^^^^-a^-{v,w)) 

^TA{b){FA{w,r2)-FA{w,0))+TB{b){FB{v,ri)-FB{v,Q)) _ ^ 

This part of the cross-section corresponds to the case when the two observed partons 
interact with each other. Obviously this requires the two partons to be produced at the 
same point in the transverse space {b = c). As a result, this part is smaller than I^^i 
corresponding to independent production part , by ~ l/A^^^ . It is important to recall 
that the independent production part has been found without expansion in powers of 
1/A, so that I^^^ fully represents the difference from independent production. 

To make the integrals over ri and convergent we change unity in the last bracket 
to exp(—T4 (6)^4(1/;, 0) — TB{b)FB{v, 0)) since this does not produce terms which could 
contribute at p, g and present the resulting cross-section in the form 

lTBiP:yuP,y2:q) = Pa{v)Pb{w) J d%TA{h)TB{h) j (iVirfV2e^^'^i+^^'-^ 

^TAmFA{w,r2)-FA{w,G)) _ ^-TA{b)FA{w,0)\ f ^mb){F(,{v,ri)-FB{v,0)) _ ^~TB{b)FB{vfl) 



j^^-TA{h)FA{w,0) I ^n(b){Fi{v,n)-FB{vfl)) _ ^-TB{b)FB{v,0) \ ^ ^-TB{b)FB{v,0) 



^TA{b){FA{w,r2)-FA{w,0)) _ ^-TA{b)FA{w,0)\ _^ ^-TA{b)FA{w,0)-TB{b)FB{v,0) | = ^ I^^B (36) 

^ ^ i=l 

The most non-trivial correlations are given by the first term. The simplest correlation 
is expressed by the last term. Indeed, dropping terms proportional to S^{p) or S^{q), 



J^) = {27:yp^{v)P^{w)d'ip + q) [ d'bTA{b)TB{b)e-^^^'^''^^^''^-^'^^'^^-^^^ 

dp J 

(37) 

and shows back-to-back correlations. The second and third terms are expressed via 
partonic and nuclear single inclusive cross-sections: 

7^1 = i2^rPB{w)^^^^^ I d%TB{b)e-^^^'^^-^-''hA{p + q, v, b) (38) 
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where we defined the single cross-section at fixed b as the integrand in (20) (for the 
central collision, /3 = 0) 



Ia{p, V, b) = PA{v)TA{b) j d\e'P'- i^e'^BmF,M-F,iv,o)) _ g-T,(fe)F,(.,o) j (39) 
The third term corresponds to v,p ^ w,q 

Jif = {2irrPAv)^^^^^ J rf^6r^(6)e-^-(^)^-(-°)/^(p + g,^,6) (40) 
Finally in terms of Ia,(b){p + Q, v{w), b) 

lT^ = j d%(^l dH^^^^^lA{p + l,v,b)lB{q-l,w,b) 

-a'''\v,w)lAip,v,b)lB{q,w,b)] (41) 



IV. CORRELATIONS 

We restrict ourselves to the study of the forward-backward multiplicity correlations. 
Since the correlations are in any case small we may restrict ourselves to a linear de- 
pendence of the average multiplicity in the backward window at a fixed multiplicity in 
the forward window function of the latter; 



{nB)nF 

{ub) 



np) = a + bnp (42) 



Here {ub) is the overall average of the multiplicity in the backward window (at all np). 
Thus defined coefficient b shows the relative deviation of the conditional average {nB)np 
when the number of jets observed in the forward rapidity window np changes by unity. 
It can be expressed via averages of linear and bilinear products of multiplicities: 

^_ 1 {nBnp) - {nB){np)^ ^^^^ 



{iib) {np) - {npY 



One gets this expression by multiplying (jl^ first by plnp) then by npplnp), summing 
over np and solving the arising system of linear equations for a and b. 

To compute the mentioned bilinear products we first point out that n! in the denom- 
inator in the phase space volume for n identical particles should refer only to particles 
within the same phase space volume. For our problem it implies that identical particles 
produced in different rapidity windows can be considered as different (see Appendix 
for details). This allows to immediately obtain simple expressions for the squares of 
multiplicities coming from independent pair production described by the product of 
single inclusive cross-sections. Consider first emission of two particles into the forward 
rapidity window. At a fixed overall impact parameter /3 we find 

{nF) = cJ ^lA{y,p) (44) 
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and 



{npinF -l)) = c j ^^^^p lAA{l3,yi,Pi,y2,P2) (45) 

where integrations over y are restricted to the forward rapidity window and c = 
l/{(faAB{P)/d'^P) where a ab{.I3) / d'^ (3 is the AB inelastic cross-section at fixed impact 
parameter j3. The latter is practically unity for (3 < Ra + Rb- As shown in section ITTl 
(Eq. (dni)), in the first approximation, the double inclusive cross-section in (45) is just 
a product of two single inclusive ones. This gives a relation 



{npinp - 1)) = ^{np)^ (46) 



which, with c ~ 1 leads to 



{nl) - {npf = {np) (47) 

(effectively the distribution seems to be Poissonian). Note that the dispersion is dif- 
ferent from zero. This is the reason why for the production of the pair into the same 
(forward) rapidity window we can limit to the first approximation in powers of l/A or 
1/B. 

Passing to the emission of jets into different rapidity windows, we may consider the 
jets different. So instead of ()45|) we shall find 

{nFnB) = cj j^-^ lAB{p,yi,p,y2,q) (48) 

From independent production (part I^^) we shall get just the product of average 
multiplicities, so that 

{ufUb) - {nF){nB) = J ^^i^|^^^4B(/3, 1/2, g) (49) 

where we used c ~ 1. 

This gives for the correlation coefficient b 

^_ J dyid^pdy2d^ql'^l{(3,yi,p,y2,q) ^^^^ 



J dyid^plAiyup) J dy2d'^qlA{y2, q) 

where integrals over yi,p and y2,q go over the forward and backward rapidity windows, 
respectively. We see that in the numerator of the expression for the correlation coef- 
ficient the leading terms in powers of 1/A and 1/B cancel and only subleading terms 
of the relative order XjA^I'^ or l/B"^/^ remain. This means that in any case forward- 
backward multiplicity correlations at a fixed (3 have the order 1/A^/^or 1/5^/^. So for 
their observation collision of comparatively light nuclei is preferable. To calculate the 
correlations one has to evaluate the integrals in ()5U|). 

Our definition of the correlation coefficient differs from the conventional one which 
reads 60 = Cov(nB, np)/ \J D{nB)D{nF) and for equal dispersions of forward and back- 
ward multiplicities is 60 = {{nFnB) — {nF){nB))/{{ny) — {nF)'^). As one can see compar- 
ing 60 to (jlH)) the difference is by factor (n|) — {nB){nB)- As this factor is proportional 
to the backward rapidity window width Ay2 (see fj47jl). the correlation coefficient P3| 
has an advantage compared to bo having no explicit dependence on the chosen rapidity 
intervals. However the value of correlation coefficient b is not limited from above by 
unity as for bo. 
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V. NUMERICAL RESULTS 

In the general case calculation of presents a formidable numerical task. To 
simplify it we limit ourselves to central collisions of identical nuclei {A = B and (3 = 0). 
We also assume the two rapidity windows narrow in rapidity and magnitudes of the 
transverse momenta |p| and |g| centered around yi,p and 1/2, q, so that in (jSUj) the 
inclusive cross-sections can be taken out of the integrals in 1/12 and \p\ and \q\ at these 
points. As to the azimuthal angle (/) between the jet momenta p and q it may be chosen 
differently in the experimental setup. If one takes into account pairs of jets with an 
arbitrary then after integration over the angles we will get a correlation coefficient 
depending only on the chosen yi,p and 1/2, 

,/ X 1 / d(f)I^^l{(3 = 0, yi,p, y2, q, 0) 

b{yi,p, y2, q) = -—, ^ 51 

71" lA{yi,P)lB{y2,q) 

We restricted our calculations to the case yi = y2 and p = q. 

(2) 

Note that the elementary processes taken into account in formally include the 
contribution from a single hard rescattering accompanied by multiple soft rescatterings. 
Physically it corresponds to smearing of the lowest order back-to-back correlation by 
soft interactions both before annd afterwards. In our formalism soft interactions are 
suppressed by a cutoff in the elementary parton-parton cross-section. However one 
can see that at large jet transverse momenta the total correlations become dominated 
by the above mentioned contribution with soft rescatterings taken at the transverred 
momenta of the order of the cutoff. Obviously such contribution cannot be desribed 
correctly in our perturbative formalism. To eleminate the effect of soft-smeared back- 
to-back correlations, especially pronounced at high momenta, we choose to restrict the 
experimental range of azimuthal angles between the jets, excluding from it angles 
close to the back-to-back configuration. We introduce a "veto angle" for the final 
momenta, that is we demand that the azimuthal angle should not be larger than 
TT — 0„eto/2. This excludes the undesirable interval of around n of length 0„eto- The 
value 71 — (pyetof^ then serves as the upper limit of integration in the numerator of (j3T|) 
and the denominator is modified by a factor of (vr — 0i,eto/2)/7r. 



For the partonic cross-sections we have taken the same expression as in [ll[ . Namely 
we assumed the effective partonic distributions with the gluon-gluon cross-section as 
a dynamical input. The lowest order cross-section was multiplied by the so-called K 
factor to take into account higher order corrections. The infrared regularization was 
realized by a non-zero gluon mass pq. The scale in the partonic distributions was taken 
as Pq + p1 where pi was the transferred transverse momentum. 

We studied S-S and Pb-Pb colisions at two cm. energies 200 and 6000 GeV. In 
accordance with [l^ we set Po = 1 GeV and i^' = 1.04 at 200 GeV and Po = 2 GeV and 
K = 2 ai 6000 GeV. In our calculations we take (pveto = 7r/6. We also present results 
for 200 GeV -A Cu-Cu and Au-Au collisions and for 200 GeV- A S-S with 0^eto = tt/S. 

In figures 1 and 2 we show the correlation coefficient 6 as a function of p at ?/i = 
7/2 = 1,2 and 3 for S-S and Pb-Pb collisions at cm. energy 200 GeV. Figures 3 and 4 
show the same correlation coefficients at 200 GeV for Cu and Au nuclei and figures 5 
and 6 represent the correlation coefficient for S and Pb at 6000 GeV. 

The found correlation coefficients drop with atomic numbers and energy and grow 
with jet momenta. For Su-Su collisions at 200 GeV they have a sharp rise for momenta 
around 7-j-8 GeV/c and rapidity ys = Vf = 3. One of the reasons for this was already 
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FIG. 1: Correlation coefficient for sulphur-sulphur collisions at s/s = 200 GeV. 
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FIG. 2: Correlation coefficient for lead-lead collisions at y/s = 200 GeV. 



discussed in this section, another is due to to kinematics: the single inclusive cross- 
sections in the denominator of (43) tend to zero at the limits of the phase space volume. 
As expected, for Pb-Pb collisions the correlation coefficients are an order of magnitude 
smaller than for S-S collisions. 
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FIG. 3: Correlation coefficient for copper-copper coUisions at s/s = 200 GeV. 
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FIG. 4: Correlation coefficient for gold-gold collisions at ^/s = 200 GeV. 



VI. CONCLUSIONS 



Our calculations show that in spite of the fact that in heavy nuclei long range mul- 
tiplicity correlations between jets are small, of the order A^"^^^ for identical nuclei, 
they are of observable magnitude if nuclei are not too heavy. They also visibly grow 
with the transverse momentum and for S-S collisions at 200 GeV. For transverse mo- 
menta of jets of about 10 -r- 12 GeV/c they stay large (10 -v- 20%) even if the back-to 
back correlations are totally excluded by 7r/3 veto angle. Their observation and mea- 
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FIG. 5: Correlation coefficient for sulphur, -s/s = 6000 GeV. 
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FIG. 6: Correlation coefficient for lead, ^/s = 6000 GeV. 
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surement for non-zero veto angles will favor a hypothesis that multiple hard collisions 
indeed occur before the fragmentation of jets into hadrons and are described by the 

perturbative QCD mechanism. However at the supposed LHC energies of 6 TeV these 
correlations are strongly suppressed reaching 0.5% for S-S collisions at jet momenta of 
about 15 GeV/c. 

Of course we understand that other effects may somewhat change our predictions 
based on the simple Glauber approach. Among them the most prominent is quenching 
of jets as they propagate through the nuclear medium. Also if the jet energy hap- 
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pens to be insufficiently fiigfi tfie formation lengtli may sfiorten to allow for the jets 
hadronization inside the nucleus, which will spoil our simple picture of hard rescatter- 
ing. Another important point is the jet fragmentation which in general can lead to 
correlation pattern for hadrons different form that for partons. These questions are 
under our consideration at present. In any case these other effects will produce certain 
corrections to the results presented in this paper, which thus may serve as a natural 
starting point. 
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VIII. APPENDIX. IDENTICAL PARTICLES IN DIFFERENT PARTS OF 
THE PHASE VOLUME (RELATED TO EQ. (47) 

Consider two particles characterized, say, by their rapidities only iji and y2 which 
take values in the interval [0,1]. The integrals over intermediate state of two identical 
particles, between which no distinction is made, in general have the form: 

i-l^'^f(m.m) (62) 

where f{y 1,1/2) is some function symmetric in yi and y2. Now we split the single 
particle phase volume in two, say [0,l] = [0,l/2] + [l/2,l] and assume that particle with 
its rapidity in [0,1/2] is particle 1 and that with its rapidity in [1/2,1] is particle 2. 
Now instead of a single intermediate state we have 3 different ones: with two particles 
1 (state A) , with two particles 2 (state B) and with one particle 1 and one 2 (state C) . 
The same integral will now be given by a sum 

Ia + Ib + Ic^ f ^^of^V (yi,y2) + / ^^^f^V (yi,y2) + / dyi f dy2f(yi,y2) 

Jo ^' Jl/2 ^' Jo Jl/2 

(53) 

As we see this sum is equal to / exactly. So considering particles in [0,1/2] and [1/2,1] 
as different gives the correct result. 
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